First Passage Properties of the Erdos-Renyi Random Graph 
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We study the mean time for a random walk to traverse between two arbitrary sites of the Erdos- 
Renyi random graph. We develop an effective medium approximation that predicts that the mean 
hrst-passage time between pairs of nodes, as well as all moments of this first-passage time, are 
insensitive to the fraction p of occupied links. This prediction qualitatively agrees with numerical 
simulations away from the percolation threshold. Near the percolation threshold, the statistically 
meaningful quantity is the mean transit rate, namely, the inverse of the hrst-passage time. This 
rate varies non-monotonically with p near the percolation transition. Much of this behavior can be 
understood by simple heuristic arguments. 
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I. INTRODUCTION 

In this article, we study a basic first-passage charac- 
teristic of random walks on Erdos-Renyi (ER) random 
graphs p|, namely, the mean time for a random walk to 
traverse between two arbitrary sites on the graph. The 
ER random graph is constructed by taking N sites and 
introducing a bond between each pair of sites with prob- 
ability p. When p = 1, all possible links exist and this 
construction gives the complete graph, where each site is 
connected to all the other N — 1 sites in the graph. As 
p decreases, the random graph undergoes a percolation 
transition at p = p c = 1/N that shares many com- 

mon features with percolation on regular lattices. An- 
other geometrical feature that is relevant for our study 
of first-passage characteristics is a second connectivity 
transition at p\ = In N/N. For p > Pi, all nodes belong 
to a single component (in the limit N — > oo), while for 
p < Pi disjoint clusters can exist 0,0- 

Much effort has been devoted to determining basic 
properties of random walks in disordered or heteroge- 
neous environments, [E IS 0i sucn as * ne ER ran dom 
graph, as well as small-world [l| and scale-free networks 
|l0| . On the random graph, there has been considerable 
work in determining how basic time scales of the random 
walk depend on the size of the graph. These include the 
mixing time - the time scale that determines how the 
probability distribution approaches its limiting behavior 
[E ITTl IT2I , the cover time - the time for a random walk 
to visit all sites of the graph ^3j; and the first-passage 
time - the time for a random walk to traverse between 
two specified points for the first time [Tpj or to return to 
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its starting point for the first time 01 ■ 

The goal of this work is to understand by simple phys- 
ical arguments how a basic first-passage property of ran- 
dom walks on random graphs depends on the concentra- 
tion of bonds in the graph. We will begin by studying 
the mean first-passage, or transit, time T a b for a ran- 
dom walk to reach an arbitrary site b on the graph when 
starting from another arbitrary site a. In general, this 
transit time T a b does not necessarily equal Tb a and it 
is sometimes convenient to consider the mean commute 
time K a b = T a b + Xf, a to avoid the asymmetry in T a b. 
However, in averaging over all pairs of sites a and b to 
obtain statistically meaningful quantities, the asymmetry 
is eliminated and (K a b) — 2(T a {,). 

In practice, mean transit and commute times diverge 
when the graph consists of more than one component 
(cluster), because a random walk that starts in one com- 
ponent cannot access sites in different components. We 
are thus led to study the dependence of the mean com- 
mute rate lZ a b = as a function of p on the random 
graph. This rate equals zero for two sites in different com- 
ponents, so that its configurational average is meaningful. 
Focusing on the commute rate is analogous to consider- 
ing mean conductance of a random conductor-insulator 
mixture near the percolation threshold, rather than the 
mean resistance. The conductance is well-behaved near 
the percolation threshold, while the resistance is diver- 
gent for all p in any finite-size system. 

In the next section, we construct an effective medium 
approximation that predicts that the mean commute 
time is independent of p for p > p±. We also find that 
the transit time varies weakly with p in the small-dilution 
limit. In Sec. Ill, we present simulation results for the 
mean commute time and the mean commute rate and 
find an unexpected non-monotonic behavior for the lat- 
ter quantity as a function of p when p w p c . In Sec. IV, 
we outline the relation between the commute time on a 
network and the two-point conductance on the same net- 
work when each link is a unit resistor. We use this electri- 
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cal network connection to explain the non-monotonicity 
of the commute rate, first for a tree structure that is a 
subset of the random graph (Sec. V), and then for the 
random graph itself (Sec. VI). 



II. ANALYTIC APPROACHES FOR THE 
TRANSIT TIME 

A. Effective Medium Approach 

We now develop an effective medium approximation for 
the commute time of a discrete-time random walk on the 
random graph [15| . In a single time step, a walk located 
at a site that is connected to z other sites can hop with 
probability 1/z to any of these neighbors. To compute 
the mean time for such a random walk to go between two 
arbitrary sites on any graph by a sequence of nearest- 
neighbor hops, we use the underlying backward equation 
[la. fl7| . This equation relates the transit time from site 
a to site b to the transit times from the neighboring sites 
of a to site b as follows: 



Tab = Pn tu = J2 Pa ^ 6t + TU 



(1) 



The first sum is over all paths II from a to b, Pn is the 
probability for the random walk to take the path II, and 
tjj is the transit time from a to & along this path. For 
each path, we then decompose the full transit time into 
the time to go from a to an intermediate site i after one 
step plus the time to go from i to b. Thus p a ^% = l/z a is 
the probability of hopping from a to i in a single step, z a 
is the degree of a, and St is the time for each step of the 
random walk. Without loss of generality, we take St = 1. 

Let us now construct an effective-medium approxima- 
tion for the average transit time on the random graph, 
under the assumption that the graph is connected. This 
condition implicitly restricts the validity of our approach 
to the range p > pi, where all nodes belong to a sin- 
gle component. A schematic representation of a random 
graph, to illustrate our approach, is shown in Fig. ^ Be- 
tween two sites a and b on the graph, a direct link to 
b may exist (thick line) with probability p. If there is 
no such direct link, then an indirect path must be fol- 
lowed. After a single step on this indirect path (medium 
lines), there may be a direct link to b with probability p 
(dashed) , or no direct link with probability q = 1 — p. 

Let us denote by r the mean transit time to go from 
a to b under the assumption that a direct link exists, 
and t' the transit time from a to b in the absence of a 
direct link. Then from Eq. Q and following an effective- 
medium assumption, r obeys the recursion formula 
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(N - l)p 
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(N - l)p 
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(2) 

The first term accounts for the walk that goes directly 
from a to b. This contribution corresponds, in Eq. Q , to 




FIG. f: Schematic decomposition of a random graph with 
starting site a and target site b. The direct link is shown as 
the thick solid line. After one step via an indirect path (to 
the sites in of the ovals), either 6 can be reached directly with 
probability p (dashed line), while with probability q = I — p 
there is still no direct connection to b. 



the case where the intermediate site i coincides with b. 
Since (TV — l)p links emanate from a on average, then ac- 
cording to the effective-medium approximation, the prob- 
ability that a random walk steps along the direct connec- 
tion is just 1 / (N — l)p. The second set of terms accounts 
for those walks in which the first step goes to an inter- 
mediate site i rather than hitting b directly. In this case, 
we again apply an effective-medium approximation and 
posit that after one step of the walk, a direct connec- 
tion from i to b exists with probability p, or no direct 
connection exists with probability 1 — p (Fig. ^) . 

To close this equation, we need an expression for r', the 
first-passage time in the absence of a direct connection 
to b. Applying the same effective-medium approximation 
as that used in Eq. J2J), we assume that after the first 
step of the walk, the terminal site b is directly reachable 
with probability p, while b is not directly reachable with 
probability q. Thus r' obeys 



r' =p(l + T)+q(l + T'). 



(3) 



Solving Eqs. @ & © gives T = N-~ and r' = N. 
Finally, we average the transit time over all pairs of ter- 
minal points and over all graph configurations. Again 
in the spirit of an effective medium approximation, this 
average is simply 



(T) = (T ab ) = pr + (1 - p) T ' = N-1. 



(4) 



Surprisingly, (T) is independent of p. Thus according 
to the effective medium approach, the complete graph 
solution, (T) = N - 1, holds for all p. 

The backward equation for the mean transit time can 
be extended to any positive integer moment of the transit 
time. Consider, for example, the mean-square transit 
time. As in the case of the mean time, the governing 
equation can formally be written as 



rpl 

1 ab — 
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Pn t 



n L n> 



(5) 



3 



For each path, we follow Eq. Q and again write the 
transit time t n as 1 + t n? , namely, the sum of the time 
for the first step and the time for the remainder of the 
path. Thus 



n 

= ^P n (l + 2*n< +t 2 n ,) 1 
n 



(6) 



In going from the second to the last line of this equation, 
we use the fact that Pn = ^2sPiPa'> where p t is the 
probability of hopping from the starting point to one of 
its nearest neighbors i, and Pn' is the probability for the 
remainder of the path II' from i to b. In the last line, the 
quantities Ti and rf are the mean and mean-square times 
to reach b when starting from i and the sum is over all 
neighbors i of the starting point. Strictly speaking, we 
should write (r) and (r 2 ) for these moments, so that it is 
obvious that (r 2 ) ^ (r) 2 . In the following, we drop these 
angle brackets because the linear and quadratic powers of 
time always appear separately and there is no ambiguity 
about where the angle brackets should appear. 

The last line of Eq. JHJ is now a backward equation for 
the second moment of the first-passage time, in which 
the previously-determined first moment is an input to 
this equation. This construction for the mean-square 
transit time can be generalized straightforwardly, albeit 
tediously, to any positive integer moment of the first- 
passage time. For the random graph, the recursion for- 
mula for the mean-square transit time is, in close analogy 
with Eqs. © and ©. 



function obeys the recursion formula 
1 N- 2 



F(z) = 



N - 1 



N- 1 



zF(z). 



This equation encodes the fact that after a single step 
(the factor z) the walk hits the target site with proba- 
bility 1/(JV-1), while with probability (N - 2)/{N - 1) 
the walk hits another interior site of the graph, at which 
point the first-passage process is renewed. The solution 
to this equation is 



F(z) 



N - 1 



1 



N 



from which 



F(t) = 



N - 1 



N-2 



N - 1 \N-1 



(7) 



Now consider the random graph with bond occupation 
probability p. Let T(t) be the first-passage probability 
from a to b when a bond is present between these two 
sites, T'(t) the first-passage probability when this bond 
is absent, and let T{£) and T'{z) be the respective gen- 
erating functions. In the spirit of our effective medium 
approximation given in Eqs. J5J and 0, we now have 
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(N~l)p 



1- 



1 



(N-l)p. 



\pzT(z) + qzT'{z)\ 



T'{z) = pzT{z) + qzT'(z). 



From these two equations, the average first-passage prob- 
ability (F{z)) — p!F(z)+qJ r '(z) has the same form as the 
first-passage probability for the complete graph (Eq. J7J)). 
Hence all moments of the transit time are independent 
of p in the effective medium approximation. 
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1 



(N — l)p I (N- l)p 

[p{l + 2t + t 2 ) + q(l + 2t' + t' 2 )] , 
t' 2 = p(l + 2t + t 2 ) + q(l + 2t' + t' 2 ). 

Using our previously-derived results for the first mo- 
ments, t — N — - and r' = N, these recursion formulae 
are easily solved. We then compute the configuration av- 
eraged mean-square transit time, (T 2 ) = pr 2 + qr' 2 , and 
obtain (T 2 ) = {2N - 3)(N - 1). Thus again, the sec- 
ond moment is independent of p and equals the second 
moment of the transit time on the complete graph. 

More generally, we show that the first-passage proba- 
bility between any two sites on a random graph, and thus 
all moments of the first-passage time, are independent of 
p in the effective-medium approximation. As a prelim- 
inary, we first compute the first-passage probability on 
the complete graph. Let F(t) be the probability that a 
random walk hits the target site for the first time at time 
t, and let F(z) = ^F(i)z* be the corresponding gener- 
ating function. For the complete graph, the generating 



B. Small Dilution Limit 

We may understand the exact dependence of the mean 
transit time in the limit p — > 1 by considering configura- 
tions with a single missing bond. There are four distinct 
cases to consider: (i) missing link between a and b (1 
configuration), (ii) missing link between a and an interior 
point (N — 2 configurations), (iii) missing link between 
an interior point and b (N ~ 2 configurations), and (iv) 
missing link between two interior points (all remaining 
configurations). Let us denote the mean transit times 
from a to b for these 4 configurations by t\, t 2 , t 3 , and 
ti. By considering each case separately, we obtain the 
recursion formulae (cf. Eq. |JT|)): 



(l+t 4 ) 

n iV-3, 
+ + j^y(1+*3) 

N - 4 

(l+*2)+jy3y(l + t4), 



h = 


l + t a 






*2 = 


1 


N- 
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with solution 

h 



N + 1 t 2 = N-l 
N t 4 = N -I 



3 

N 
2 

N' 



(8) 



Then, by averaging over the appropriate number of con- 
figurations for each class, we obtain the mean transit time 



(T) 



N 3 - 2N 2 + N + 4 
N(N — 1) 



N- 1 + 



4 



(9) 



To interpret this result, note that the absence of a 
single bond corresponds to a bond concentration p « 
1 — 2/N 2 . At this value of p, the mean transit time has 
the asymptotic behavior (T)/(T(p=l)) - (1 + 4/7V 3 ) ~ 
[1 + V2(l ~ p) 3/2 ]- The first correction to (T) is thus 
of the order of (1 — p) 3 ^ 2 , rather than linear in (1 — p), 
as one might naively expect. This small first correction 
to (T) near p = 1 makes plausible the effective-medium 
result that (T) is independent of p. 




FIG. 2: Mean commute time (dashed) and mean commute 
rate (solid) versus bond occupation probability p for a random 
graph of N = 100 sites. Both quantities are normalized to 
have the value 1 for the complete graph (p = 1). Averages 
over 10 3 graph realizations were performed for each p. Also 
shown are the locations of p c — 0.01 and pi w 0.046. 



III. SIMULATION RESULTS 

To test the effective-medium prediction for the mean 
first-passage time, we now turn to numerical simulations. 
For very small systems (N < 8) , we have obtained the 
exact first-passage time by averaging over all configura- 
tions of random graphs, over all pairs of endpoints, and 
over all random walks. For the graph configuration av- 
erage, each realization is weighted by the factor p k q E ~ k , 
where k is the number of occupied links in the graph, 
E = N(N — l)/2 is the total number of possible links, 
and q = 1 — p. We then average over all pairs of end- 
points directly. By this averaging, the mean transit time 
is simply one-half of the mean commute time. Rather 
than averaging over individual walks directly, we solve 
exactly the recursion formulae in Eq. for the transit 
times between all pairs of points. 

For larger systems, the exact enumeration of all graph 
configurations is impractical. Instead we average over a 
finite number of graph realizations and endpoint pairs, 
but still performed the exact average over all random 
walk trajectories by numerically solving Eq. (Q. For ef- 
ficiency, we start our simulation with an empty graph, 
add bonds one at a time and then update the commute 
times between all pairs of sites in the graph after each 
bond addition. Each graph is then weighted by p k q E ~ k 
so that we can obtain the commute time as a function 
of p. We repeat this sequential graph construction over 
many realizations. The graphs that we obtain by this 
sequential growth are the same as those obtained by a 
static construction in which each bond is present with 
probability p = 2M /N(N - 1) when N is large (see Q 
for different, but equivalent ways of constructing random 
graphs). 



For the average commute time, we only include con- 
nected graphs in the ensemble, while for the average rate, 
the ensemble consists of all graph configurations. This re- 
striction plays a significant role only for p < p\, where 
the random graph normally consists of multiple compo- 
nents. Typical results for a graph of 100 sites are shown 
in Fig. [5] Above the connectivity threshold p\ — In N/N, 
the average transit time varies slowly with p, in agree- 
ment with our effective medium approach. The apparent 
singularity of the average commute time at a value p < p\ 
stems from finite size effects. 
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FIG. 3: Mean commute rate Kona random graph for N — 50 
(o), N = 100 (A), 200 (□), 400 (o), and 800 (v) sites as a 
function of the average site degree fi = p(N — 1). These rates 
are normalized to one for the complete graph limit. Averages 
over 10 3 graphs were performed for each case. 
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The behavior of the mean commute rate is shown in 
Figs. [21 and [3J Unexpectedly, this rate is non-monotonic 
in p for p w p c = 1/N, as shown in detail in Fig. For 
this plot, we use the average degree, [i — p(N — 1) as the 
dependent variable, because it has the desirable feature 
that the percolation transition occurs at the same value 
fi c = 1 for all N. The fact that the non-monotonicity in 
the commute rate occurs near /j, = 1 suggests that this 
anomaly is connected with the percolation transition of 
the random graph. 

To understand this non-monotonicity, we first make a 
connection between the commute rate and the conduc- 
tance on the same network when each link is a unit re- 
sistance, and then analyze the structure of the random 
graph in the critical regime to evaluate the conductance. 
This result will then be used to infer the dependence of 
the commute rate on the mean degree \x. 

IV. ELECTRICAL NETWORK CONNECTION 



10 2 F 
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FIG. 4: Average two-point conductance (thick solid curve) 
on a random graph with N = 100 sites. The dotted line 
corresponds to G = /Lt/2, the asymptotic large- p form for the 
conductance. 



In principle, first-passage properties of random walks 
on a graph can be obtained from the underlying Lapla- 
cian matrix of the graph 5]. The eigenvalue spectrum 
of the Laplacian provides many time-dependent random 
walk characteristics. This matrix formulation also re- 
veals deep analogies between random walks on a graph 
and the electrical network problem on the same graph in 
which each occupied link is a resistor of unit resistance 
(see Ref. IS] for a nice exposition of these connections). 

For example, the commute time K a b between a and b 
and the conductance G a b between these same two sites 
are simply related by K a b = 2M/G a b @,IH- Here M is 
the number of bonds in the cluster that contains both 
a and b. Equivalcntly, the mean commute rate lZ a b = 
I /Kab is given by 



Ttnh — 



Gab 

2M' 



(10) 



As we shall see, it is much easier to estimate the conduc- 
tance rather than the commute rate of a random graph 
by direct means. We will then rely on this connection 
between G a b and lZ a b to determine the latter quantity. 

For reasons of numerical convenience, we will often 
consider the following sum of the rates 



Tl a 



(11) 



We include the factor of 2 in the definition because TZ a 
then equals 1 for the complete graph. We may also sum 
freely over all sites b in the system Eq. (|ll(l because 
G a b = for any sites that are not in the same cluster 
as a. Finally, we obtain the average commute rate for 
the graph by averaging over all initial sites a: 



(12) 



In the limit of large fj, all the sites in the graph belong 
to the same cluster and M = fiN/2. Thus the average 
commute rate becomes 



ah 



(13) 



where G is the two-point conductance averaged over all 
pairs of graph endpoints. Thus, as we have discussed, 
first-passage times and two-point conductances are inti- 
mately connected. 

For the conductance itself, it is worth noting that this 
function behaves anomalously near the connectivity tran- 
sition. Although the conductance must increase mono- 
tonically with /x |l§j , the rate of increase changes for /i in 
the critical range between 1 and In AT (Fig. HJ. For large 
[/,, the conductance asymptotically approaches G = fx/2 
(dashed line), a result that corresponds to the average 
commute rate approaching 1Z = 1, in agreement with the 
result of Fig. El 



V. STRUCTURE OF THE RANDOM GRAPH 

To determine the behavior of the commute rate for 
general values of fi, we first need to resolve the structure 
of random graphs at these values of fi. From this struc- 
tural information, it is relatively easy to determine the 
conductance, from which we may then infer the behavior 
of the commute rate. 

An advantage of formulating the average commute rate 
as in Eq. Ijlll) is that only sites that are connected to the 
starting point a contribute to lZ a . Thus we can restrict 
the endpoint b to lie in the cluster that also contains a. 
Then averaging over many realizations of these clusters 
is equivalent to averaging over a in Eq. (|12|l . 
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A. Rooted Geodesic Tree 

To generate a cluster within the random graph that 
contains the starting site a of the random walk, we first 
construct a subset of the cluster that we term the rooted 
geodesic tree (RGT). We can then build the rest of a ran- 
dom graph cluster from the RGT. The RGT is a specific 
subset of a random graph cluster that: (i) spans all the 
sites in the cluster, and (ii) the distance between a and 
any site b on the RGT is also the shortest distance be- 
tween these two sites in the random graph cluster. The 
notion of the RGT is inspired, in part, by the minimal 
spanning tree, a construction with useful applications in 
network flow problems 19, 20]. 




FIG. 5: A rooted geodesic tree (RGT). Disks represent sites. 
The largest black disk is the root. Disks in successive shells 
have successively smaller radii. The solid lines are links in the 
RGT. The broken lines are bonds that are added subsequently 
to generate a random graph cluster. 

To construct the RGT, we start with TV sites and no 
bonds, and assign one site (denoted by a) as the root; 
this site is defined to be at level j = 0. We generate the 
RGT as a series of successive shells centered about a. The 
j th shell, denoted by Sj, contains those sites that can be 
reached from a after exactly j hops between connected 
neighbors. We define Sj as the expected number of sites 
in the j th shell. 

Suppose that we have just generated the j th shell. An 
unassigned site y becomes part of shell Sj+i if a link is 
created that joins y to an arbitrary site x in Sj. For each 
unassigned site, at least one such link will be created with 
probability 1 — (1 — p) Sj . We then test each possible link 
between y and the sites in Sj one by one. Each such link 
will be created with probability p. When the first such 
link is created, y becomes an element of Sj + i and we then 
consider the next unassigned site for potential inclusion 
in Sj+i. If none of the possible links is created, y re- 
mains available for inclusion in subsequent shells. After 
all the unassigned sites have been tested, the shell Sj+i is 
complete. This growth process continues until either no 
unassigned sites remain or if all attempts to incorporate 
the available sites into the current shell fail. In the latter 
case, the total number of sites in the RGT is less than 



N. 

There are two important subtleties associated with this 
construction algorithm for an RGT. First, bonds that are 
not examined in the initial construction of the RGT can 
only exist between sites in the same or in adjacent shells 
of the RGT. A second important point is that in building 
the RGT, each examined bond was tested one time only 
and is therefore included in the RGT with probability p. 

Since the number of sites in successive shells of the 
RGT grows exponentially in the number of steps away 
from the root, the radius of the largest cluster (giant 
component), is given by the criterion [i L w (3N. Here 
(3 — /3(ju) is the fraction of the initial N sites that belongs 
to the largest cluster. By the definition of the RGT, the 
fraction of sites in the largest cluster in the random graph 
and in the underlying RGT are identical. Thus the radius 
of the RGT is given by 



L 



In AT 
In \x 



ln/3 
In \x 



(14) 



Since (3 is a rapidly growing function of /i 0, 0, the 
radius of the giant component of the RGT is an increasing 
function of fj, just above the percolation threshold ji c = 1. 
This increase in radius occurs because the RGT acquires 
progressively more sites with increasing fi. On the other 
hand for sufficiently large /i, the giant component will 
contain almost all sites in the graph and the radius of 
this component will decrease as \x is increased still further 
(Fig.[f)J|. This non-monotonic behavior of the RGT radius 
on [i is ultimately connected to the non-monotonicity in 
the commute rate. 




FIG. 6: Average RGT radius (solid curve) as a function of the 
average degree fi. The dotted curve is the average commute 
rate on the same structure; this rate has been scaled to lie on 
a similar vertical range. Data are based on 10 5 realizations of 
RGTs of 100 sites. 



B. Role of Loops on Graph Structure 

Given a rooted geodesic tree, it is possible to augment 
the tree to generate a realization of a random graph clus- 
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ter. We merely attempt to add to the RGT each of the 
bonds between sites on the cluster that were not previ- 
ously considered in the construction of the RGT itself. 
Each such bond addition attempt is carried out with 
probability p. As a result of the fact that each of these 
newly-added bonds and each bond in the RGT is present 
with probability p, all bonds in the full random graph 
cluster are present with probability p. Furthermore, both 
the RGT and the corresponding random graph have the 
same number of sites and radius for the same values of 
fi and N. As a result of this equivalence, the RGT un- 
dergoes the same percolation transition as the random 
graph itself when fi passes through 1. 

VI. THE COMMUTE RATE 

We now use the connection between commute rate and 
conductance to understand the non-monotonicity in the 
commute rate for a random walk on a random graph. 
As indicated in Fig. there are three regimes for the 
commute rate: (I) an initial increase with /i for small fi] 
(II) a decrease over an intermediate range; and (III) an 
ultimate increase for large /i. For regimes I and II, the 
commute rates on the RGT and the random graph are 
nearly identical and it is simpler to consider the commute 
rate on the RGT. We then investigate how adding the 
links to the RGT to create a random graph affects the 
commute rate. 



A. Commute Rate on the RGT 

For a tree graph, the resistance between two sites is 
simply the path length between these two sites. Thus 
the average commute rate in Eq. has the form 

b^a j = l 

Here the number of links in a tree is one less than the 
total number of sites V, and T> a b is the distance between 
a and b. Thus TZ a is the inverse moment of the distance 
between the root a and all other sites in the tree. The sec- 
ond equality follows from the shell structure of the RGT, 
where L is the radius of the tree. Thus we need only the 
statistics of the shell sizes of the RGT to determine the 
commute rate. 

Since each realization of the RGT is distinct, the num- 
ber of sites V, the radius L, and the shell sizes Sj fluc- 
tuate from realization to realization. To calculate the 
configuration-averaged commute rate (1Z a ), we first use 
the algorithm of the previous section to generate RGTs. 
Then we solve the random walk problem on each realiza- 
tion and average Eq. I|15|) over realizations to determine 
the commute rate (Fig. [7J . 

To understand the non-monotonicity of the commute 
rate for the RGT, consider first the small-/j, limit. Bc- 




0.5 1 1.5 2 2.5 3 3.5 4 

FIG. 7: Commute rates on the RGT (o) and on random 
graphs (A) for N = 100 sites based on averages over 10 5 
realizations. The lower curve is our RGT prediction lZ a = 
ln/z/lniV for the large-/x limit. The upper curve in the pre- 
diction for the random graph TZ = (/j, — (Eq. 1171 '). The 
approximate locations of regimes I, II, and III are indicated. 



cause isolated sites contribute zero to the average rate, 
the commute rate must initially increase with /i, as small 
trees begin to form. Once most sites are no longer iso- 
lated, the radii of typical RGTs then increase with [i due 
to the merging of small trees. This increase in radius 
causes a decrease in the commute rate, as can be seen by 
writing the rate in Eq. (|15f) as 

n a {L) = Et/x> ( 16 ) 

3=1 3 I 3 = 1 

In the analogous expression for lZ a {L + 1), the numer- 
ator increases by Sl+i/(L + 1) while the denominator 
increases by Sl+i- Thus lZ a (L) is a decreasing function 
of L, so that a tree with a larger radius will have a smaller 
commute rate. 

As argued in Sec. V, a further increase in /i will cause 
the radius of the RGT to eventually decrease with jU. 
Correspondingly, the commute rate enters regime III and 
increases with fi. In this regime, we now use the fact that 
the number of sites in successive shells of the RGT grows 
exponentially in the distance from the root. Thus the 
shell at radius L contains almost all of the sites of the 
RGT. As a naive approximation, we then replace the sum 
in Eq. (|15fl by the last term to give, in the limit of large 

(Ka)RGT ~VL h[N- 

This result agrees extremely well with numerical results 
for the commute rate on RGTs, as shown in Fig. [7] 
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B. Role of Loops on Commute Rate 

We now investigate how adding loops to the RGT to 
build a random graph affects the behavior of the mean 
commute rate. Starting with a realization of an RGT we 
generate a cluster of the random graph by adding missing 
bonds, following the procedure discussed in Sec. V. The 
addition of these bonds will create loops that provide 
alternative paths between the root site and the endpoints 
of a random walk (Fig. |SJ. The ostensible effect of these 
additional paths is to increase the commute rate between 
the root and any endpoint. 



b 




FIG. 8: Schematic representation of the random graph. The 
included RGT is also shown. Loops typically arise at a dis- 
tance L — \nN/lnfi from the root. A site in this last shell 
will typically have fi independent paths to the root. 



from the root are at the same potential. Thus the con- 
ductance of the two joining RGTs is simply one-half of 
the conductance between the root and the last shell of a 
single RGT. 

For this last step, we approximate the RGT by an in- 
finite Cayley tree with branching ratio fi. The resistance 
between the fc-th and the (k + l) st shell in this tree is 
H~( k+1 ' , since the links between the two shells are in 
parallel. Because the shells are in series, the resistance 
from the center to infinity is simply the geometric sum, 
Sfc°=o^ = Tjzp Thus the conductance between 
a and b is G a t, = (/i — l)/2. Substituting this result in 
Eq. {EH) ; then gives the commute rate 



This result converges to 1 as /i — > oo, in agreement with 
the effective medium approach in Sec II as well as our 
simulation results. Closer to the percolation threshold, 
however, Eq. i|17|) and simulation results quantitatively 
disagree because our naive picture for the structure of 
the random graph no longer applies. 

Thus we observe that the eventual increase in the com- 
mute rate (regime III) stems from the combined effect of 
the decrease in the radius of the underlying RGT embed- 
ded within a random graph cluster and the emergence of 
loops that join two RGTs in the random graph. 

VII. CONCLUSION 



To estimate the two-point conductance for a random 
graph for general \i > In N, we start with the picture that 
the graph consists of two RGTs, one emanating from a 
and the other from b (Fig. EJ). For a graph of N sites, 
the radius of each tree is of order L ~ ln(iV/2)/ln/x. We 
argue that these two trees tend to join only at the out- 
ermost shell because this is where most of the sites in 
the trees are located. We further assume that, in the 
equivalent resistor network, all sites at the same distance 




FIG. 9: Schematic random graph structure to calculate the 
conductivity between two sites a and b for mean degree fj, > 1. 
An RGT is grown around both a and b. The two RGTs meet 
at a distance 0(ln N/ In fj,) from each of a and b. Broken 
lines are links between sites in the outermost shells of the two 
respective RGTs. 



We studied a basic first-passage characteristic of ran- 
dom walks on random graphs that is related to the time 
for a walk to travel between two arbitrary points on a 
graph. We first constructed an effective medium theory 
and a small dilution approximation for this mean transit 
time. The former approach predicted that the mean tran- 
sit time, and also all positive integer moments of the tran- 
sit time, are independent of the bond concentration p for 
p greater than the connectivity threshold p\ = hi N/N. 
The small dilution approximation also predicts a slow 
dependence of the transit time on p near p = 1. Our nu- 
merical simulation results are in qualitative accord with 
a transit time that is slowly varying in p for p > p\ . 

Below the connectivity threshold, the transit time is 
not well defined because the mean time for a random 
walk to hop between sites on different components of 
a disconnected graph is infinite. To avoid this pathol- 
ogy, we studied the inverse of the commute time, namely, 
the commute rate. We developed a simple heuristic pic- 
ture for the behavior of the commute rate that relied on 
first identifying an embedded rooted geodesic tree (RGT) 
within an arbitrary random graph cluster. For the RGT, 
it is simple to compute the commute rate in terms of a ge- 
ometric picture for the tree and thus argue that this rate 
is a non-monotonic function of p in the critical regime. 

We then presented a simple physical picture for the 
influence of loops on the behavior of the commute rate. 
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Qualitatively, the dependence of the radius of the un- 
derlying RGT on the bond concentration explains the 
behavior of the commute rate close to the percolation 
threshold. For larger fi loops become an important fac- 
tor and are ultimately responsible for the non-monotonic 
dependence of the commute rate on p. While our argu- 
ments were heuristic and the approximations made are 
uncontrolled, they provide an intuitive picture for the 
structure of random graphs and also provide qualitative 
and satisfying agreement with simulation results for the 
commute rate. 



Acknowledgments 



We thank Paul Krapivsky and Federico Vazquez for 
helpful suggestions. VS and SR gratefully acknowledge 
financial support from NSF grants DMR0227670 (at BU) 
and DOE grant W-7405-ENG-36 (at LANL). DbA sim- 
ilarly acknowledges NSF grant PHY0140094 (DbA) for 
financial support. 



[1] P. Erdos and A. Renyi, Publ. Math. Debrecen 6, 290 
(1959). 

[2] B. Bollobas, Modern Graph Theory (Springer, New York, 
1998). 

[3] S. Janson, T. Luczak, and A. Ruciriski, Random Graphs, 
(John Wiley & Sons, Inc., New York, 2000). 

[4] B. D. Hughes, Random Walks and Random Environ- 
ments. Volume 2: Random Environments, (Clarendon 
Press, Oxford, 1996). 

[5] L. Lovasz, in Combinatorics, Paul Erdos is Eighty, Vol. 
2 (eds. D. Miklos, V. T. Sos, T. Zsonyi), Janos Bolyai 
Mathematical Society, Budapest, 1996, 353-398. 

[6] W. Woess, Random Walks on Infinite Graphs and 
Groups, (Cambridge University Press, Cambridge, 2000). 

[7] D. Aldous and J. A. Fill, Reversible Markov 
Chains and Random Walks on Graphs. http://stat- 
www .berkeley.edu/users/aldous/ RWG /book. html. 

[8] D. ben-Avraham and S. Havlin, Diffusion and Reactions 
in Fractals and Random Systems, (Cambridge University 
Press, Cambridge, 2000). 

[9] E. Almaas, R. V. Kulkarni, and D. Stroud, Phys. Rev. E 
68, 056105 (2003). 
[10] J. D. Noh and H. Rieger, Phys. Rev. Lett. 92, 118701 
(2004). 



[11] A. J. Bray and G. J. Rodgers, Phys. Rev. B38, 11461 
(1988). 

[12] R. Monasson Eur. Phys. J. B 12, 555 (1999). 

[13] C. Cooper, A. Frieze, Proceeding of the fourteenth annual 

ACM- SIAM symposium on Discrete algorithms, 140- 

147, (SIAM, 2003) 
[14] N. Masuda and N. Konno, Phys. Rev. E, 69, 066113 

(2004). 

[15] For general presentations of effective medium theory, see 
e.g., S. Kirkpatrick, Rev. Mod. Phys. 45 , 574 (1973); 
J. Koplik, J. Phys. C 14, 4821 (1981). 

[16] N. G. van Kampen, Stochastic Processes in Physics and 
Chemistry, 2nd ed. (North-Holland, Amsterdam, 1997). 

[17] S. Redner, A Guide to First-Passage Processes (Cam- 
bridge University Press, New York, 2001). 

[18] P. G. Doyle and J. L Snell, Random Walks and Elec- 
tric Networks. Caruso Mathematical Monographs #22 
(Mathematical Association of America, Oberlin, OH, 
1984); also math.PR/0001057. 

[19] See R. K. Ahuja, T. L. Magnanti and J. B. Orlin, Network 
Flows: Theory, Algorithms and Applications, (Prentice 
Hall, Englewood Cliffs NJ, 1993). 

[20] See also Sec. 5.2 of Ref. Q for a construction related to 
the RGT. 



